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Abstract. We review the relation of Af=4 superconformal multi-particle models on the real 
line to the WDVV equation and an associated linear equation for two prepotentials, F and U. 
The supcrspacc treatment gives another variant of the integrability problem, which we also 
reformulate as a search for closed flat Yang-Mills connections. Three- and four-particle 
solutions are presented. The covector ansatz turns the WDVV equation into an algebraic 
condition, for which we give a formulation in terms of partial isometrics. Three ideas for 
classifying WDVV solutions are developed: ortho-polytopes, hypcrgraphs, and matroids. 
Various examples and counterexamples are displayed. 
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1 Introduction 



x 



Over the past decade, there has been substantial progress in the construction of AA=4 supercon- 
formal multi-particle mechanics (in one space dimension) [ll[2j[3j|4l[5l[6l[7l[8l[9l [10] . In 2004 
a deep connection between these physical systems and the so-called WDVV equation [UJ H] 
was discovered [4j, relating Calogero-type models with D (2,1; at) superconformal symmetry to 
a branch of mathematics concerned with solving this equation [13|, [T41 [To] [To] [T71 [T~8] 119]. Here, 
we describe physicists' attempts to take advantage of the mathematical literature on this subject 
and to develop it further towards constructing and classifying such multi-particle models. 

There exist different versions of the WDVV equation in the literature, so let us be more 
specific. Originally, the WDVV equation appeared as a consistency relation in topological field 
theory, where the puncture operator singles out one of the coordinates, so that the associa- 
ted Frobenius algebra is unital and carries a constant metric |11[ fT2[ I13j . A few years later, 
a more general form of the WDVV equation appeared as a condition on the prepotential F 
of Seiberg-Witten theory (four-dimensional j\f=2 super Yang-Mills theory) |20[ I21 [ [22]. Here, 
the distinguished coordinate is absent, and so the Frobenius structure constants do not lead to 
a natural metric. This so-called generalized WDVV equation takes the form 



A j A k 



Vi, j,k = 1,... ,n 



(1) 



*This paper is a contribution to the Proceedings of the Workshop "Supersymmetric Quantum Me- 
chanics and Spectral Design" (July 18-30, 2010, Benasque, Spain). The full collection is available at 
|http://ww.emis.de/joumals/SIGMA/SUSYQM2010.html] 
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for a collection of n x n matrix functions (Aj)^ m = — didid m F. However, any invertible linear 
combination of these matrices yields an admissible metric 

rj = rf Ai so that A i r]~ 1 Aj = AjT]^ 1 A i . (2) 

i 

It is easy to see that this metric can be absorbed in a redefinition [14} 115]. 
Ai := r]~ 1 Ai — > [A i ,A j ]=Q with ^rf A = 1, 

giving a formulation equivalent to (PQ). For a constant metric, e.g. r/ = A\ = const, we fall 
back to the more special case which arose in topological field theory. Since 1999, Veselov and 
collaborators have been constructing particular solutions to ([2]), introducing so-called V-systems 
and featuring a constant metric r] = £\ x 1 A{ [151 El CEH1 [19] . 

In comparison, AA=4 supersymmetric multi-particle models are determined by two prepo- 
tentials, F and U, the first of which is subject to the generalized WDVV equation [3]. Here, 
the conformal invariance imposes a supplementary condition on our matrices, which amounts to 
choosing the Euclidean metric 

^ ^ x Aj = 1 y Ai = Ai, 

i 

and so one may drop the label 'generalized'. The map to Veselov's formulation is achieved by 
a linear coordinate change, x l — > x^M^ with r] = MM ± [23] , 

The goal of the paper is fourfold. First, we would like to review the appearance of the 
WDVV equation in the construction of one-dimensional multi-particle models with sti(l,l|2) 
symmetry [H EJ EJ 13 El [10] - In particular, we draw the attention of the mathematical readers 
to the second prepotential U, which enlarges the WDVV structure in a canonical fashion, and to 
the superspace formulation, which yields an alternative formulation of the integrability condition. 
Second, we plan to provide some explicit three- and four-particle examples for the physical model 
builders. Third, we intend to rewrite Veselov's notion of V-systems in a manner we hope is more 
accessible to physicists, using the notion of partial isometry and providing further examples. 
Fourth, we want to advertize some novel attempts to attack the classification problem for the 
WDVV equation. The standard ansatz for the propotential employs a collection of covectors, 
which are subject to intricate algebraic conditions. These relations may be visualized in terms 
of certain polytopes, or in terms of hypergraphs, or by a particular kind of matroid. Neither of 
these concepts is fully satisfactory; the classification problem remains open. However, for low 
dimension and a small number of covectors they can solve the problem. 



2 Conformal quantum mechanics: Calogero system 

As a warm-up, we introduce n + 1 identical particles with unit mass, moving on the real line, 
with coordinates x l and momenta pi, where i = 1, 2, . . . , n + 1, and define their dynamics by the 
Hamiltoniar0 

H = y iP i + VB(x\...,x n+1 ). (3) 
For the quantum theory, we impose the canonical commutation relations (H = 1) 
[x\pj] = i<5/. 



1 Equivalently, it describes a single particle moving in R" +1 under the influence of the external potential Vb- 
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Together with the dilatation and conformal boost generators 

D = — \{x l pi +PiX l ) and K = \x % x % , 
the Hamiltonian (J3J) spans the conformal algebra so(2, 1) in 1 + dimensions, 

[D, H] = -iH, [H, K] = 2iD, [D, K] = iK, 

if and only if {x l di + 2)Vb = 0, i.e. the potential is homogeneous of degree —2. If one further 
demands permutation and translation invariance and allows only two-body forces, one ends up 
with the Calogero model, 

v B = y 92 

B ^ (x i - xi) 2 
i<j v > 

3 A/"=4 superconformal extension: su(l, 1| 2) algebra 

Our goal is to A/"=4 supersymmetrize conformal multi-particle mechanics. The most general 
A/"=4 extension of so(2, 1) is the superalgebra D(2, 1; a), but here we specialize to D(2, 1; 0) ~ 
sw(l, 1|2) 3 su(2). Further, we break the outer su(2) to by allowing for a central charge C. 
The set of generators then gets extended 



(H, D, K) -> (H,D,K,Q a ,S a ,J a ,C) with a = 1,2 and a = 1,2,3 

and hermiticity properties (Q Q )^ = Q a an d 05a) = S a . 
The nonvanishing (anti)commutators of stt(l,l|2) read 

[D, H] = -iH, [H, K] = 2iD, 

[D, K] = +iK, [J a , J b ] = ie abc J c , 

{Q a , Q 13 } = 2H5J, {Q a , S?} = +2i(a a )Jj a - 2DSJ - iCSj, 

{S a ,S?} =2K5j, {Q a ,S f} } = -2i(a a )p a J a -2D6p a + iC5p a , 

[D, Q a ] = —-^Qai \D,S a ] = +2*Sa> 

[K, Q a ] = MS a , [H, S a ] = -iQa, 

[Ja,Qa] = -\{aa) a ^Qp, [J a ,S a ] = -\{<J a) a ^ S p , 

[D,Q a ] = -^Q a , [D,S a ] = +^S a , 

[K,Q a ] = +iS a , [H,S a ] = -iQ a , 

[Ja, Q°] = W(°a)8 a i [Ja, S°] = ^(a a )g a . 



To realize this algebra on the (n+l)-particle state space, we must enlarge the latter by adding 
Grassmann-odd degrees of freedom, tp a and ^) %a = ip a \ with i = 1, . . . , n + 1 and a = 1, 2, and 
subject them to canonical anticommutation relations, 

{Vi, 4} = 0, {i> ia ,i> j ?} = and {C V^} = Sjfi. 
In the absence of a potential (subscript '0'), the generators are given by the bilinears 

Qo a =Pi^l Qo=Pd ia and S 0a = x i ^ S% = x^ ia , 

H = \piPi, D = -l^pi+pix*), K = \x i x\ J 0a = \^ iot {a a )J^, 
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where a a denote the Pauli matrices. Surprisingly however, the free generators fail to obey the 
su{\,\\2) algebra, and interactions are mandatory! The minimal deformation touches only the 
supercharge and the Hamiltonian, 

Qa = Qoa-i[Soa,V], Q a = Q% - i[S Q , V] and H = H + V, 

keeping 5 = So, S = So, D = Do, K = Kq and J = Jo. 

Being a Grassmann-even function of ip, ip and x, the potential V may be expanded in even 
powers of the fermionic variables. It turns out that we must go to fourth order for closing the 
algebra, i.e. [HUE] 

V = V B (x) - UyixMW ) + \F\ m (x)Wrf a P%), (4) 

where the angle brackets (• • • ) denote symmetric (or Weyl) ordering. The functions Uij and Fijki 
are totally symmetric in their indices and homogeneous of degree —2 in {x 1 , . . . , x n+1 }. For 
completeness, we also give the interacting supercharge, 

Q a = ( Pj - ix%(x))^ - ±x l F ijkl (x){^ k ^ l a ). (5) 

4 The structure equations for (F, U): 
WDVV, Killing, inhomogeneity 

Inserting the minimal ansatz @ for V into the sit(l,l|2) algebra and demanding closure, one 
finds that 

Uij = didjU and F ijM = didjd k diF 

are determined by two scalar prepotentials U and F, which are subject to so-called structure 
equations [21 HJ E] , 

(didkdpF^dpd^F) - (dAdpF^dpdkdjF) = o, x%d j d k F = s jk , (6) 

didjU - (didjd k F)d k U = 0, x i d i U = -C. (7) 

The left equations (|6h,) and ([Th.) are homogeneous quadratic in F (known as the WDVV equa- 
tion) [1 1|, I12j and homogeneous linear in U (a type of Killing equation). The right equations ([Hb) 
and (|7b) introduce well-defined inhomogeneities, so that the prepotential must be of the form 

F = - \x 2 In x + F hom and U = -C In x + U hom (8) 
with Fhom of degree —2 and Uh om of degree in x. This also shows the redundancies 

U — U + const and F ~ F + quadratic polynomial, 
which for F is also apparent in the twice-integrated form of ([6b), 

(x%-2)F = -ixV. 

It is convenient to separate the center-of-mass dynamics from the relative particle motion, 
since the two decouple in all equations. The center-of-mass motion is already nonlinear but 
explicitly solved by ([8]) without homogeneous terms (the central charge is additive). In new 
relative- motion coordinates, which again we name x % but with i = 1, 2, . . . , n, the configuration 
space is reduced to M n . The Killing-type equation ([Th) implies, as its compatibility condition, 
the WDVV equation ([6^l) contracted with djU. Furthermore, the contraction of ([6h.) with x l is 
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trivially valid, thanks to ((6)3) . This effectively projects the WDVV equation to n — 1 dimensions. 
Since its symmetry is that of the Riemann tensor, it comprises as many independent equations, 
namely ^^(n — l) 2 (n — 2) in number. In particular, ([6^) is empty for up to three particles and 
a single condition for four particles. 

The leading part of the potential is also determined by U and 

V b = l(diU)(diU) + ^ (didjd^F) (pidjd k F) > 0, 

and the expressions in ([5]) simplify to 

x l F ijM = -djd k diF and x l U i:j = -djU. 

Therefore, finding a pair (F,U) amounts to defining an su(l,l|2) invariant (n + l)-particle 
model. For more than three particles, however, this is a difficult task, and very little is known 
about the space of solutions. 

5 Superspace approach: inertial coordinates in M. n+1 

When analyzing super symmetric systems, it is often a good idea to employ superspace methods. 
This is also possible for the case at hand, where the construction of a classical Lagrangian seems 
straightforward in M=A superspace [25] 126] l27l l28| [29] . 

For each particle, we introduce a standard untwisted N= A superfield 

m A (t,6 a ,6 a ) = u A (t) + 0(6,6) with A = 1, . . . , n + 1, 

obeying the constraint^] 

D 2 u A = = D 2 m A -> d t [D a ,D a ]m A = — > [D a ,D a ]m A = 2g A 

with constants g , which will turn out to be the coupling parameters. The general jV=4 super- 
conformal action for these fields takes the form@ 

5 = - J &td 2 6d 2 6G{vL) = \ J dt [G A b{u)u a u b - G A B{u)g A g B + fermions] 

already written in [I], with a superpotential G(u) subject to the conformal invariance condition 
G - G A u A = \c A u A 

for arbitrary constants c A , so that it is of the form G = — ^culnu + terms of degree one. 

Generically, such sigma-model-type actions do not admit a multi-particle interpretation, how- 
ever, unless the target space is flat. This requirement imposes a nontrivial condition on the 
target-space metric Gab(u) [S], 

Riemann(G j4 s) = < — ► G A[BX G XY G YC]D = 0. 

Equivalently, there must exist so-called inertial coordinates x l , with i = 1,2,... ,n + l, such that 

S = J dt \bijX l x? — Vq(x) + fermions . 
2 Here and later, we sometimes reinstate h to ease the interpretation. 

3 The constants g A can be SU(2)-rotated into the constraints, so that D 2 u A = ig a = -D u A but [D a , "D a ]u A = 0. 
Subscripts on G denote derivatives with respect to u, i.e. G A = dG/du A etc. 
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The goal is, therefore, to find admissible functions u = u (x) and compute the corresponding 
G and Vq. The above flatness requirement leads to a specific integrability condition for u A := 



diU , namely 

o^aH x )) = {( u °»)~ 1 ) a ='■ W ^,i = 9 i™A = -q-t(x), (9) 

which says that the transpose of the inverse Jacobian for u — > x is again a Jacobian for a map 
w — > x. This defines a set of functions wa{x) dual to u a (x), in the sense that their Jacobians 
are inverses [9], 

WA,iU B i = S A B < > W A ,iU A j = Sij. 

Equivalent versions of the integrability condition ([9]) are [9] 

u[A i d o uB } = < — > W[A,i9jWB],i = 0, (10) 
fijk '■= —WA,idkU A j is totally symmetric, 
fijk = didjdkF and fi m [kfl]mj = 0) 

which includes the WDVV equation for F. In contrast, there is no formulation purely in terms 
of U. 

Conformal invariance restricts u A to be homogeneous quadratic in x, hence wa to be homoge- 
neous of degree zero (including logarithms!), thus f^k is of degree —1. The second prepotential U 
is also determined by u A (x) via 

U(x) = —g A WA{x) so that C = —x l diU = CAg A , 

and automatically fulfills the Killing- type equation (J7^). The classical bosonic potential then 
reads 

y B cl = \(d$J){diV) = \g A g B w A , l w B , l . 
Finally, for the superpotential G(u) the integrability condition becomes 

u^u^Gab = Sij < — > Gab = -w A ,iW B ,i = -d A w B = -d B w A , 
so that, up to an irrelevant u- linear shift of G, 

wa = —G a < — > G = —u a wa, 
and we have 

U = g A G A and h 3 k = -\u A u B u c k G A BC -> V£ l = -\g A g B G AB ■ 

However, knowing the superpotential does not suffice: the relation between x 1 and u A is needed 
to determine U(x) and F(x). On the other hand, if a solution F to the WDVV equation can be 
found, this problem reduces to a linear one [9]: 

u % + hjku k = or w A,ij ~ fijkWA,k = 0, with f ijk = didjd k F. (11) 

Finally, we remark again that the center-of-mass degree of freedom can be decoupled, so that 
all indices may run form 1 to n only. 
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6 Structural similarity to closed flat Yang— Mills connections 

It is instructive to rewrite our integrability problem in terms of n x ra-matrix-valued differential 
forms, in a compact formulation closer to Yang-Mills theory. To this end, we define 

{vfy := u, {-didjF) := f and (-didjd k Fdx k ) := A = A k dx k . (12) 
Since A k = d k f and didjd k F = WA,idkU A j, we have 

A = df ->■ dA = and A = u^du -)• dA + AAA = 0, 
from which we learn that 

= AAA = ±d[/,d/] = -du- 1 Adit, (13) 

which is nothing but the WDVV equation again. Hence, we are looking for connections A 
which are at the same time closed and flat. Dealing with a topologically trivial configuration 
space M n , it implies that A is simultaneously exact and pure gauge. The exactness is already 
part of the definition (|12p . and the pure-gauge property is what relates A with u. We remark 
that A and / are symmetric matrices while u is not. Furthermore, the inhomogeneity ([6b) 
demands that x l dif = 1. The task is to solve (fl~3j) for / and for it, which then yield d 3 F and 
VU = -2u~ x g. 

Of course, we cannot 'solve' the WDVV equation by formal manipulations. But even given 
a solution A (and hence /), it is nontrivial to construct an associated matrix function u. For 
this, we must integrate the linear matrix differential equation (jlip . 



dw T = A« T , (14) 
which qualifies u as covariantly constant in the WDVV background. The formal solution reads 

oo 

u T = ^2 f (k) with /(0) = i, fW=f an d d//< fc > = d/C fc+1 \ 

k=0 

up to right multiplication with a constant matrix. The matrix functions are local because 
d(d// (fc) ) = -df A d/ (fc) = -df A d//^ 1 ) = -A A A/^" 1 ) = 

due to the WDVV equation. Likewise, one has 
/d/( fc) = /d//^ 1 )=d(//( fc )-/( fc+1 )). 

Note that the naive guess u T = is wrong since [/,d/] = d(/ 2 — 2f^) ^ 0. 
We provide two explicit examples for n = 2, with the notation 

i=l i=2 j 2,2 2 

x =: x, x =: y and x + y =: r . 
Starting from the B2 solution with a radial term [6j [9] 

F = — \x 2 lnx — ^y 2 lny — \{x + y) 2 ln(x + y) — \ {x — y) 2 ln(x — y) + \r 2 lnr, (15) 
we have 

l(H(x 2 -y 2 )£] In ^{ 

2 I In ^ l n [( x 2_ y2) ^ 
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with (xd x + yd y )f = 1 and, hence, 

A df ^ X<2 ~ ^ 2 [ ^ X ^ | + 4x 2 y 2 f x dx — ydy xdy — ydx\ 
xyr 4 \ xdy J (x 2 — y 2 )r 4 \xdy — ydx xdx — ydy I 

It is easy to check that indeed A A A = but [A, f] ^ 0. The solution to (|14p turns out to be 



U = — T 

r 4 



/ xr 4 


yr 4 \ 


r=i 


\xy 4 


yx A j 





u 1 = \r 2 , 

u 2 = \x 2 y 2 /r 2 



with an arbitrary non-degenerate constant matrix T, as may be checked by inserting it into (|14p . 
One may also begin with a purely radial WDVV solution [9], 

2 2 

zji i 2i . s i n 2\-,i , x ~~ y , x y 

F = -^r\nr — ► / = _(lnr )1 + <r 3 + ^-<7i, 



2x 2y \ r=i j u 1 = r 2 , 



and find 

u = r , 

2:r arctan f - ?/ 2y arctan f + lu 2 = r 2 arctan f . 

For more generic weight factors in (|15p , u 2 is expressed in terms of hypergeometric functions [9] . 



u 



7 Three- and four-particle solutions 

An alternative method for constructing solutions (F, U) attempts to find functions u A (x) sat- 
isfying (fTOl) . It is successful for n + 1 = 3 since the WDVV equation is empty in this case. 
Imposing also permutation invariance, a natural choice for three homogeneous quadratic sym- 
metric functions of (x*) = (x,y,z) is 

(x + y + z) 2 , 
i 2 = {x-yf + {y- z) 2 + (z-x) 2 , 

u 3 = [(2x -y- z)(2y -z- x)(2z -x- y)] 2/3 h(s), (16) 

where h is an (almost) arbitrary function of the ratio 

[(2x — y — z)(2y — z — x)(2z — x — y)] 2 
[(x - y) 2 + (y - z) 2 + (z — x) 2 ] 3 

Not surprisingly, (|16p fulfils the integrability condition ()10p . so we are guaranteed to produce 
solutions. It is straightforward to compute the Jacobians u\ and WA,i and proceed to the 
prepotentials. Writing (g A ) = (51,52,53), the bosonic potential comes out as 



T/ cl 



5?/24 



(x + y + z) 2 
1 



+ 



1 

324 



(1 - 2s)g 2 + 2s 



(h + 3sh') 2 



(x - y) 2 

(x + y + z) 2 
5(2 + 35) 



+ 



+ 



(y - z) 2 (z - x) 



9l 



4sl- s g 2 g 3 + 2s l 3~ 25 g 2 



324(1 + 



+ 



9l 



[l + 3J) 2 (x - yf + (y- z) 2 + {z 



+ 



+ 



(x - y) 2 (y - z) 2 (z - x) 



W=4 Multi-Particle Mechanics, WDVV Equation and Roots 



9 



where in the second equality we specialized to 
h(s) = s 5 <— > 



5 3 _ [(2x -y-z)(2y-z-x)(2z-x- y)] 2 ^ +2S 



[(x - y) 2 + (y- z) 2 + (z- x) 2 ] 3S 
Putting 53 = for simplicity, the corresponding prepotentials are 

U = ~ ln(x + y + z) - 18 ^ g ^_ 3§ j H x ~y)(y- z){z - x) 

-^^Hix-yf + iy-zf + iz-xf], 

F = -\(x + y + z) 2 ln(s + y + z) - \ [(x - y) 2 ln(x -y) + (y- zf hx(y - z) 
+ {z- x) 2 \n{z - x)] + [(2a -y-z) 2 ln(2x - y - z) 
+ (2y — z — x) 2 ln(2y — z — x) + (2z — x — y) 2 ln(2z — x — y)] 

+ ! |> - v? + (y- + (« - x ?\ MO* - y? + {y- + - *) 2 ] • 

We recognize the roots of G2 plus a radial term in the coordinate differences. The potential 
simplifies in two special cases: 

5 = 44> h = 1 : ^bHs'i = 53 = 0) is P ure Calogero, 

S = I O h = s 1 / 6 : Vn{g\ = 52 = 0) is pure Calogero. 

In the full quantum potential, Vq = + tL-F 1 " F" 1 , the couplings g" 4 receive quantum correc- 
tions. 

Stepping up to four particles, i.e. n + 1 = 4, it becomes much more difficult to construct 
solutions, since the integrability condition is no longer trivial. Our attempts to take a known 
WDVV solution and exploit the linear equations (jlip for u \ have met with success only spo- 
radically. In most cases, the hyper geometric function 2 -Pi turns up in the expressions. A simple 
permutation-symmetric example uses the A% solution with a radial term, 

for which we discovered [9] 
it 1 = (x + y + z + w) 2 , 

u 2 = (x - y) 2 + (x- z) 2 + {x- w) 2 + {y- z) 2 + {y- w) 2 + (z - w) 2 , 
u 3 = u 2 I ( x + y-*- w \ and .4 = U 2 I (P_ 



with 

P 2 



(x - y + z - w) + 2\J (w - x)(y - z), q 2 = (x - y - z + w) + 2-sJ (w - y)(x - z) 
At 



and I(x) = 

Jo 



The Jacobians and the bosonic potential are algebraic but not of Calogero type. It remains 
a challenge to find (m 2 ,u 3 ,u 4 ) for the ^3 WDVV solution without radial term, 



^ i ' i i<j 
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8 Covector ansatz for prepotential F 

For the rest of the presentation, we concentrate on the WDVV equation in M n , 

(did k d p F)(d p did 3 F) - (d i d l d p F)(dj,d k djF) = o with (x% - 2)f = -±xV, 

since, together with £7=0, its solutions already produce genuine A/"=4 superconformal mechanics 
models. Leaving aside a possible radial term 

F rad = -r 2 lnr with r 2 := ^{x*) 2 , (17) 

i 

we employ the standard 'rank-one' or 'covector' ansatz [2] 
F = — \ y~](q • x) 2 In a ■ x 

a 

containing a set {a} of covector s 

a = (ai,a2, ■ ■ ■ ,a n ) E (M n )* or S i(M n )* — > a{x) = a ■ x = onx 1 , 
subject to the normalization 

^ajOj = Sij < — > y^ct (gi a = 1 (18) 

a a 

which takes care of the inhomogeneity in (|17p . The WDVV equation turns into an algebraic 
condition on the set of covectors [HI [151 E] > 

^ — — — (aiipj - ajPi)(a k {3i- atiPk) = with a • /3 = 5 lj aif3j. (19) 

Apart from the normalization (|18p . the covectors are projective, so we may think of them as a 
bunch of rays. Let us denote their number (the cardinality of {a}) by p. We may assume that 
no two covectors are collinear. Since an orthogonal pair of covectors does not contribute to the 
double sum, two mutually orthogonal subsets of covectors decouple in (fT9|) . and it suffices to 
consider indecomposable covector sets. In n = 2 dimensions, ([T8]) implies (fl~9j) . but already for 
the lowest nontrivial dimension n = 3 only partial results are known (2[ IT?] [T8| [T9l [6] . 



9 Partial isometry formulation of WDVV 

Let us gain a geometric understanding of (|19|) . Each of the \p(p — 1) pairs (a,/3) in the double 
sum spans some plane it ~ aA/3 G A 2 ((IR n )*), but not all of these planes need be different. 
When we group the pairs according to these planea^l, the tensor structure (a A f3)® 2 of (fT9|) tells 
us that this equation must hold separately for the subset of coplanar covectors pertaining to 
each plane ir, 

y °-/?Mfl 2 =0 v.. (20) 

' a ■ x B ■ x 



5 A given covector may occur in different pairs, thus in different groups. Covectors are not grouped, only their 
pairs. 
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Depending on the number q of covectors contained in a given plane tt, one of three cases oc- 
curs ng 

case (a) tt contains zero or one covector — > equation trivial, 

case (b) tt contains two covectors, tt ~ a A f3 — > orthogonality a ■ /3 = 0, 

case (c) tt contains q > 2 covectors — )• projector condition on tt : 

^ a a = A^IU =: A^ for A T eR and P 2 = P w with rank(P 7r ) = 2. (21) 

The latter is the proper covector normalization for the planar subsystem, which implies the 
(trivial) WDVV equation on tt to hold. Establishing the projector condition (|2ip simultaneously 
for all planes is a nontrivial problem, since covectors usually lie in more than one plane, which 
imposes conditions linking the planes. 

For a more quantitative formulation, we express (|2ip in terms of partial isometries. After 
introducing a counting index a = 1, . . . ,p for the covectors {a} = {a\, . . . , a p }, we collect their 
components in an n x p matrix A. This defines a map 

A: R p ^R n given by A = (a^)^'^ with AA T = t n , 

encoding the total normalization (118|) . For each nontrivial plane tt, we select all a aB € tt, 
s = 1, . . . , q, via 

B w : W^RO by {a a } ^ {a J 
and write the combination 

A w : R» by Ar := = (n,„Jl ^ 

Our projector condition then reads 

A^Al = KP n <— ► AjA 7r =A 7r Q 7r (22) 

with projectors on M n and Q,,- on R 9 of rank two and multipliers A^, for any nontrivial 
plane tt. Therefore, A is a WDVV solution iff -4= is a rank-2 partial isometry (|22p for each 
nontrivial plane 7r! An alternative version of (|22p is 

^7T ^7T ^7r = ^TT^-tt • 

Note that A ^ A n B n . Since the projectors are of rank 2, we may split A n over R 2 : 

3 Dj, : R q -> IR 2 and CV : M 2 <- IP such that A n = CjD w . 
The situation can be visualized in the following noncommutative diagram: 
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We illustrate the partial isometry formulation with the simplest nontrivial example, which 
occurs at n = 3 and p = 6, by providing a one-parameter family of covectors {a, f3, 7, a',f3', j'}(t) 
via 



A 



1 

6 



Q 

/ 6t 




(3 
-St 
3y/3t 





7 

-3* 
-3V3t 





with iu = V2 - 3i 2 . (23) 



It is easily checked that AA T = 1. A quick analysis of linear dependence reveals that 12 of the 
15 covector pairs are grouped into 4 planes of 3 pairs each, leaving 3 pairs ungrouped. 3 coplanar 
pairs imply 3 coplanar covectors, hence there are 4 nontrivial planes containing q = 3 covectors, 
namely 

{a (3 7), (a /3 V), («W>, (a'P'l), 
and 3 planes containing just two covectors, which are indeed orthogonal, 

a ■ a' = ■ = 7 • y = 0. 
Let us test the projector condition (I22p for two of the planes: 
'2t 



.4 



(a /3 7 > 



.4 




-f 










=> kaz 





J 




3w 


-3^^;^ 




\/3u> 


y/Sw 


1 AX 


2^/3 


2VSj 






(a/3'7') 



where the matrices on the right are idempotent. Hence, in both cases, A n is proportional to 
a partial isometry, with a (parameter-dependent) multiplier A ff . The other two nontrivial planes 
work in the same way. We have proven that (|23j) produces a family of WDVV solutions. This 
scheme naturally extends to include imaginary covectors as well. 



10 Deformed root systems and polytopes 

It is known for some time [144 E] that the set $ + of positive roots of any simple Lie algebra (in 
fact, of any Coxeter system) is a good choice for the covectors. So let us take 

{a} = <]?t = <£+ U with a L • a L = 2 and a s • «S = 1 or |. 

where the subscripts 'L' and 'S' pertain to long and short roots, respectively. Having fixed the 
root lengths, we must introduce scaling factors {f a } = {/lj/s} m 



F = - 1 2\h E E V a " x)2ln 



a • x \. 



The normalization condition (|18|) has a one-parameter solution, 



./L = ^ + (^-^ V )i, 

/ L e E ) a ® a = 1 - *• { 1 with * G 

«6*£/ |/ s = -^ + (/i-r/i v )t, 
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where h and h v are the Coxeter and dual Coxeter numbers of the Lie algebra, respectively. The 
roots define a family of over-complete partitions of unity. Amazingly, all simple Lie algebra 
root systems obey (|20p . and they do so separately for the pairs of long roots, for the pairs of 
short roots and for the mixed pairs, of any plane ir. This leads to the freedom (t) to rescale 
the short versus the long roots and provides a one-parameter family of solutions to the WDVV 
equation [HI [THE]. (In the simply-laced case there is only one solution, of course.) 
For illustration we give two examples. Let {e^} be an orthonormal basis in R n+1 . For 

A n © A\ : {a} = < — ej, ej|l<i<j<n + l> we find 

i 

F An^ = J>< - x>? ln(x* - *»■) - ^(Ez'YlnfZy) 

i<j ^ i ' ^ i ^ 

with center-of-mass decoupling, while for the non-simply-laced case (n = 2, p = 6) without 
center of mass 

G 2 ■ {a} = {^( e i - e i)» 7g( e i + e j ~ 2e k) (hj, k) cyclic j one gets 

F G2 = _IZ|^( X 1 _ x 2)2 ^ ^! _ ^ _ + x 2 _ 2x 3)2 ln ^1 _ ^2 _ ^3) 

+ cyclic. 



A natural question is whether one can deform the Lie algebraic root systems by changing 
the angles between covectors but keep (|20p valid. So which deformations respect the WDVV 
equation? Based on a few examples, we conjecture that the (suitably rescaled and translated) 
covectors should form the edges of some polytope in M. n . Non-concurrent pairs of edges then 
have no reason to be coplanar with other edges, thus better be orthogonal. Concurrent edge 
pairs, on the other hand, always belong to some polytope face, hence automatically combine 
with further coplanar edges to a nontrivial plane ir. The hope is that the polytope's incidence 
relations take care of the WDVV equation, e.g. in the form of (I22p . For p > ^n(n + 1), there is 
enough scaling freedom to finally arrange the normalization (|18p with {fa}- 

This expectation is actually bourne out in the case of the A n root system, 
which, with p = \n(n + 1), is in fact the minimal irreducible system in each 
dimension n and uniquely fixes {fa}- Starting with an arbitrary bunch of 
\n{n + 1) rays in M n , we reduce the freedom in their directions by impos- 
ing firstly the n-simplex incidence relations and secondly the orthogonality 
conditions for skew edges. Let us do some counting of moduli (minus global 
translations, rotations and scaling): 






ray moduli 


incidences 


simplex moduli 


orthogonality final moduli 


# 


\n 2 {n - 1) 


-i(n-2)(n 2 - 1) 


i(n-l)(n + 2) 


-i(n-2)(n + l) n 


n = 2,3,4 


2, 9, 24 


0,-4,-15 


2, 5, 9 


0,-2,-5 2,3,4 



We find that the moduli space M(A n ) of these so-called orthocentric n-simplices is just n- 
dimensional. It can be shown [23] that indeed it fits perfectly to a family of WDVV solutions 
found earlier [171 \18\ \19\ , lending support to our polytope idea. We remark that the previous 
example (|23|) represents a one-parameter subset in ^(^3). 

Let us make this observation more explicit in the case of n = 4. Using the recursive con- 
struction of orthocentric n-simplices presented in [6] for n = 4 and computing the corresponding 
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scaling factors is feasible but algebraically involved. Therefore, we just present a 'nice' one- 
parameter subfamily of solutions, with t € M+ and w 2 = t 2 — |, 



A 





fwy/2 


f V2 




-fv 7 ^ 


fV2 




-±V2 





\ 


1 





-f\/6 fx/6 


fV5 


fv^ 


fV5 











2i 





^^3 







f v 7 ^ 












I o 














t 









For t 2 = | we have the root system of A4, at t 2 



the first six covectors disappear and 



leave A\. When < t 2 < |, the first six covectors are imaginary, and in the singular limit t 2 — >0 
we obtain the A3 roots and fundamental weights, but can no longer maintain our normalization. 

A more familiar parametrization embeds the A4 root system into M 5 , in the hyperplane 
orthogonal to the center-of-mass covector e%, with s € K + and u 2 = 20s 2 — 10s + 1, 



A 



(l-4s)v / 5 



/« 





u 





« 





1 - 


- s 


— s 


— s 




—u 








u 
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s 


1 - s 


— s 


— s 





u 


— it 








—u 




s 


— s 


1-s 


— s 





—u 





— 1i 


—u 







s 


— s 


— s 


1-s 


Vo 

















4s - 


- 1 


4s - 1 


4s - 1 


4s- 1/ 



Now s = yields the roots of A4, beyond s = ^(1 — the first six covectors turn imaginary, 

and the singular limit s — > ^ (-u — >• |) gives the A3 roots and fundamental weights, orthogonal 
also to &i — 5e§. This pattern generalizes to an interpolation between the A n roots and the 
A n _i roots and fundamental weights. 

What about deformations of other root or weight systems? We give two 
more prominent examples in n = 3 dimensions. First, consider the p = 9 
positive roots of B3 and observe that, from four copies of them, we can 
assemble the edges of a truncated cube. It is possible to deform the latter 
into a truncated cuboid while keeping the orthogonalities and producing 
a six-parameter family of covectors, 

{a • x} ={dix 1 , (I2X 2 , d^x 3 ; C3(c2X 1 ±cix 2 ), c\(c3X 2 ±C2X 3 ), C2(cix 3 ±C3X 1 )} , Ci,d{ £l. 

The normalization ^ Q f a a®a = 1 can be achieved with 




{/«} 



c 2 + c 2 



r 2 



cj + c 2 



~ 2 d\ 



z 2 d 2 



3 2 d\ 



c 2 + c 2 + c 2 + c 2 . 



c 2 c 2 ' c 2 c 2 ' c 2 c 2 . 



One sees that the relevant combinations \ff^a depend only the three ratios ^. It turns out 
that we have constructed a three-dimensional moduli space of WDVV solutions [TT l fl8| . fl9] . 

Second, again using A3, it is possible to combine four copies of its three 
positive vector weights with six copies of its four positive spinor weights 
to the edge set of a rhombic dodecahedron, with each rhombic face being 
dissected into two triangles. There exists a three-parameter family of defor- 
mations in line with the orthogonalities, given by 



a ■ x = dix 1 , 
a + /3 + 7 



(3 ■ x 
a — (3 



-- d 2 x 2 , 
7 



7 • x ■ 



d 3 x 3 ; 
1 




-a 



P + l 
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and re-normalization is achieved by 

_d| + 4 + 4 4-4 + 4 4 + 4-4 

Ia cP4 ' Ip d?4 ' h d 2 4 ' 

2 2 2 2 
/spinor = d = + <^2 ^ ^3 • 

In this case, the combinations \ffaOi depend only on the ratios -4, and we again discover a two- 
dimensional family of WDVV solutions |18|, I19j . It seems that indeed the polytope's incidence 
relations imply the WDVV equation, thus allowing us to construct solutions F purely geomet- 
rically, by guessing appropriate polytopes with certain edge multiplicities. 



11 Hyper graphs 



Sadly, our ortho-polytope concept fails, as may be seen from the first counterexample at (n,p) 
(3,10): 




4 


4 




5 

V3 



6 7 
-y/2 

V3 -1 
3 



V6 
-1 
3 

IT 




(24) 



is properly normalized, AA = I3, and may be checked 
to fulfil the partial-isometry conditions (|22D for each non- 
trivial plane. It turns out, however, that there exists no 
polyhedron whose edges are built from (suitably rescaled 
copies of) all ten column vectors in (|24[) . In the attempt 
shown to the left, one of the would-be edges (labelled '9') 
runs inside the convex hull created by the others. 
This lesson demonstrates that it may be better to restrict ourselves to the essential feature 
of A, which is the coplanarity property of its columns a. Even though there is not always an 
ortho-polytope, we may still hope that each n x p matrix A with 



AA 1 



t, 



and 



j(a a Aa b ) 



a;, = Va, b = 1, 



■P 



already obeys the crucial conditions (|22[) for all q > 2 planes. 

Suppose we have m nontrivial planes and label them by p, = 
1, . . . ,m. The > 2 covectors in the plane 7fy are grouped in the 
subset {a a ^} C {a a }, with s = 1, . . . ,q^. A shorter way of encoding 
this coplanarity information is by using only the labels rather than 
denoting the covectors. Thus, we combine the a>s for each nontrivial 
plane ir^ in the subset {a^ \s = 1, . . . , q^} =: 11^ C {1, . . . ,p}, and then 
write down the collection H(A) := {Ili,ri2, . . . ,n m } C V({1, . . . ,p}) 
of these (overlapping) subsets. Such subset collections are known as 
simple hypergraphs [30] , They are graphically represented by writing 
a vertex for each covector label and then, for each ji, by connecting all vertices whose labels 
occur in ILj. The resulting graph has p vertices and m connections 11^, called hyperedges. 
Note that a vertex represents a covector, and a hyperedge stands for a (nontrivial) plane, 
thus gaining us one dimension in drawing^. As an example, the hypergraph for (|24f) reads 




3 Our simple hypergraphs contain only q M -vertex hyperedges with g M > 2, hence no one- or two-vertex hyper- 



edges. 
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{{1234}{1580}{2670}{179}{289}{356}{378}{457}{468}{490}} and is represented above (with 
'0'='10'). To the mathematically inclined reader, we note that our simple hypergraphs are not 
of the most general kind: they are also 

• linear: the intersection of two hyperedges has at most one vertex (uniqueness of planes) 

• irreducible: the hypergraph is connected (the covector set does not decompose) 

• complete: when adding the g At = 2 planes, each vertex pair is contained in a hyperedge 

• orthogonal: a nonconnected vertex pair is 'orthogonal' (property of the q^ = 2 planes) 

Of course, two hypergraphs related by a permutation of labels are equivalent. Thus, our pro- 
gram is to construct, for a given value of p, all orthogonal complete irreducible linear simple 
hypergraphs and check the partial-isometry conditions (|22p for each plane ir. Unfortunately, 
this is not so easy, because the orthogonality is not a natural hypergraph property but depends 
on the dimension n of a possible covector realization. In fact, it is not guaranteed that such 
a realization exists at all. Therefore, the classification of complete irreducible linear simple hy- 
pergraphs with p vertices has to be amended by the construction of the corresponding covector 
sets in W 1 , subject to the orthogonality condition. 

12 Matroids 

Luckily, there is another mathematical concept which abstractly captures the linear dependence 
in a subset of a power set, namely the notion of a matroid \31\ l32"j 133] . There exist several 
equivalent definitions of a matroid, for example as the collection {C^} of all circuits C 
{1, . . . ,p}, which are the minimal dependent subsets of our ground set {1, . . . ,p}: 

• The empty set is not a circuit. 

• No circuit is contained in another circuit. 

• If C\ ^ Ci share an element e, then {C\ U C<2)\{e\ is or contains another circuit. 

Of course, we identify matroids related by permutations of the ground set. 

The idea is that each circuit corresponds to a subset of linearly dependent covectors. Indeed, 
every n x p matrix A produces a matroid. However, the converse is false: not every matroid is 
representable in some M n . If so, it is called an M-vector matroid, with rank r < n. The rank 

= \C^\ — 1 of an individual circuit is the dimension of the vector space spanned by its 
covectors. Excluding one- and two-element circuits qualifies our matroids as simple. It may 
happen that two rank-ci circuits span the same vector space, for example if they agree in d of 
their elements. Hence, it is useful to unite all rank-d circuits spanning the same d-dimensional 
subspace in a so-called d-flat F^, with 2 < d < r. We call such a d- flat minimal if it arises 
from a single circuit, i.e. |i^| = d + 1. In this way, we may label the matroid more efficiently 
by listing all 2-flats, 3-flats etc., all the way up to r — 1. Needless to say, we are only interested 
in connected matroids, i.e. those which do not decompose as a direct sum. Also, for a given 
dimension n we study only R- vector matroids of rank r = n and ignore those of smaller rank, 
since they can already be represented in a smaller vector space. Finally, we need to implement 
the orthogonality property. So let us call a matroid orthogonal, if any pair of covectors which 
does not share a 2-flat is orthogonal. Note that further orthogonalities (inside 2-flats) may be 
enforced by the representation. 

A matroid of rank r can be represented geometrically in R r_1 as follows. Mark a node for 
every element of the ground set (the covectors). Then, connect by a line all covectors in one 

2- flat, for all 2-flats. Next, draw a two-surface containing all covectors in one 3-flat, for all 

3- flats, and so on. We illustrate this method on two examples, the A4 and the B3 matroid: 
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The A4 case has r 
1 < i < j < 5. Then, 



4, and it is natural to label the ten covectors by pairs (ij), with 



{C^} = {{(ij)(ik)(jk)},{(ij)(ik)(j£)(k£)^ 

{(lOG)fi)ffl(W)}} with (j) = (y) r(ii) 

lists ten circuits of rank 2, fifteen circuits of rank 3 and twelve circuits of rank 4. The former 
represent ten 2-flats, the middle unite in triples to five 3-flats and the latter combine to the 
trivial 4-flat, 

{F 2 } = {{(mk)(jk)}}, 

{F 3 } = {{(ij)(ik)(i£)(jk)(j£)(k£)}}, 

{F 4 } = {{(12)(13)(14)(15)(23)(24)(25)(34)(35)(45)}}. 

Orthogonality is required between pairs with fully distinct labels. The B 3 example is of rank 
three but less symmetric. We label the three short roots by i and the six long ones by i and i, 
with i = 1,2,3, and obtain sixteen rank-2 circuits grouping into seven 2-flats and thirty-nine 
rank-3 circuits combining into the unique 3-flat (i 7^ j 7^ k 7^ i), 

{C^} = {{ijk}, {ijk}, {ijj}, {ijk}, {ijk}, 

{ij ij}, {ij ij}, {ij ij}, {i i ij}, {i i ij}, {i ijk}, {i ijk}, {i ijk}, {i ijj}}, 
{F 2 } = {{1,2,3,3}, {1,3,2, 2}, {2,3, 1,1}, {1, 2,3}, {1, 2,3}, {1,2,3}, {1,2,3}}, 
{F 3 } = {{1,2,3,1,2,3,1,2,3}}. 

Here, we see that i 1 i and i _L i, but the realization in R 3 actually enforces i _L j as well. 

The task then is to classify all connected simple orthogonal R-vector matroids for given 
data (n,p). There are tables in the literature which, however, do not select for orthogonality. 
Another disadvantage is the fact that matroids capture linear dependencies of covector subsets 
at any rank up to r, while the WDVV equation sees only coplanarities. Therefore, it is enough 
to write down only the 2-flats, which brings us back to the complete irreducible linear simple 
hypergraphs again. Still, the advantage of matroids over hypergraphs is that they provide 
a natural setting for the orthogonality property and the partial-isometry condition (|22p . Once 
we have constructed a parametric representation of an R-vector matroid as a family of n x 
p matrices A, we may implement the orthogonalities and directly test (122h for all nontrivial 
planes it. A good matroid is one which passes the test and thus yields a (family of) solution(s) 
to the WDVV equation. 

Another bonus is the possibility to reduce a good matroid to a smaller good one by graphical 
methods. The two fundamental operations on a matroid M are the deletion and the contraction 
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of an element a € {1, 
these look as follows 



,p} (corresponding to a covector). In the geometrical representation 



• deletion of a, denoted M\{a}: remove the node a and all minimal <i-flats it is part of 

• contraction of a, denoted M/{a}: remove the node a and identify all nodes on a line 
with a, then remove the loops and identify the multiple lines created 

Both operations reduce p by one. Deletion keeps the rank while contraction lowers it by one. 
On the matrix A, the former means removing the column a (corresponding to the covector a a ) 
while the latter in addition projects orthogonal to a a - Connectedness has to be rechecked after 
deletion, but simplicity and the M-vector property are hereditary for both actions! Furthermore, 
contraction preserves the orthogonality, but deletion may produce a non-orthogonal matroid. 
Since the contraction of a good matroid corresponds precisely to the restriction of V-systems 
introduced by [l8j|T9], we are confident that it generates another good matroid. A similar state- 
ment holds for the multiple deletion which produces a V-subsystem in the language of |18| [T9] . 

The first nontrivial dimension is n = r = 3, where simple matroids (determined by {F2}) 
are identical with complete linear simple hypergraphs (given by {H^}). Their number grows 
rapidly with the cardinality p: 



number p of covectors 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


how many simple matroids? 





1 


2 


4 


9 


23 


68 


383 


5249 


232928 


28872972 


of the above are connected 











1 


3 


12 


41 


307 


4844 


227612 


28639649 


of the above are M-vector 











1 


3 


11 


38 


7 


? 


? 


? 


of the above are orthogonal 














1 


1 


1 


1 


3 


? 


? 



Below, we list all good (/) and a few bad (j) cases up to p = 10, with graphical/geometric 
representation and the name of the corresponding root system. Parameters s, t, u indicate con- 
tinuous moduli. 



{{123}, {145}} y ^ 
{{123}, {1456}} y 
{{123}, {145}, {356}} 



A 



M{6} 

5 3 \{4,5,9} 

D(2,l;a)\{7} 




{{123}, {145}, {356}, {246}} = ^A^ 

{{123}, {145}, {356}, {347}, {257}, {167}} 

{{123}, {145}, {356}, {347}, {257}, {167}, {246}} 
{{123}, {145}, {356}, {347}, {257}, {248}, {1678}} 
{{123}, {145}, {347}, {257}, {2489}, {1678}, {3569}} 
{{150}{167}{259}{268}{456}{479}{480}{1234}{3578}{3690}} 
{{179}{289}{356}{378}{457}{468}{490}{1234}{1580}{2670}} 



M-vector but not orthogonal ^ 
M-vector but not orthogonal ^ 
M-vector but not orthogonal I 
A 3 (s,t,u) / 
604 of ^3 = D(2,l;a)(s,t) / 
Fano matroid - not M-vector ^ 
B 3 \{9}(s,t) / 
B 3 (s,t,u) / 
C AB(l,3)(t) / 
C AB(l,3)(t) / 




The last two lines (with p = 10) arise from restrictions of a one-parameter deformation of the 
p = 18 exceptional Lie superalgebra AB (1,3) root system [18]. More precisely, the first of these 
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two cases is rigid and can also be obtained from the Eq roots, while the second case retains the 
deformation parameter. 

We have developed a Mathematica program which automatically generates all hypergraphs 
subject to the simplicity, linearity, completeness and irreducibility properties up to a given p. 
Furthermore, hypergraphs that admit no orthogonal covector realization are ruled out, thereby 
drastically reducing their number. For a generated hypergraph we then gradually build a pa- 
rametrization of the most general admissible set of covectors whereby it turns out whether 
the hypergraph is representable. A major step forward would be to completely automate this 
process also; we are confident that this is feasible. Finally, on the surviving families A(s,t, . . .) 
of covector sets, the program tests the partial-isometry property (122p equivalent to the WDVV 
equation, for all nontrivial planes n. 

A natural conjecture is that our class of hypergraphs or ma- 
troids always produces WDVV solutions, rendering this final test 
obsolete. However, running the program for a while reveals a coun- 
terexample at (n,p) = (3, 10), given by the hypergraph to the right. 
In this diagram, the hollow nodes indicate additional orthogonality 
inside a plane spanned by four covectors. We must conclude that 
a geometric construction of WDVV solutions is still missing. 

Although the connected simple orthogonal M-vector matroids are not classified and the 
WDVV property does not automatically follow from such a matroid, this approach is still useful 
in exhausting all covector solutions for a small number of covectors at low dimension, i.e. for 
a limited number of particles. In this way one of us has, in fact, proven [34[ [35] that there are 
no other four-particle solutions (n = 3) with p < 10 beyond those determined in |18[ [19]. The 
matroid itself does not capture the moduli space of solutions with a given linear dependence 
structure, but its systematic realization by an iterative algorithm will do so (as it did for n = 3). 
Around a given solution, the local moduli space may be probed by investigating the zero modes 
of the WDVV equation linearized around it. 

13 Summary 

We begin by listing the main points of this article: 

• A/"=4 superconformal n-particle mechanics in d = 1 is governed by U and F 

• U and F are subject to inhomogeneity, Killing-type and WDVV conditions 

• a geometric interpretation via flat superpotentials gave new variants of the integrability 

• there is a structural similarity to flat and exact Yang-Mills connections 

• the general 3-particle system is constructed, with three couplings and one free function 

• higher-particle systems exist, tedious to construct; hypergeometric functions appear 

• the covector ansatz for F leads to partial isometry conditions with multipliers A,,- 

• finite Coxeter root systems and certain deformations thereof yield WDVV solutions 

• certain solution families admit an ortho-polytope interpretation 

• hypergraphs and matroids are suitable concepts for a classification of WDVV solutions 

• the generation of candidates can be computer programmed 

• not all connected simple orthogonal M-vector matroids are 'good' 

There remain a lot of open questions. First, can our hypergraph/matroid construction pro- 
gram detect new WDVV solutions not already in the list of [181 IS]? Second, given a 'good' 
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matroid, can we generate its moduli space, e.g. by linearizing the WDVV equation around it? 
Third, the explicit Hamiltonian of the A^=4 four-particle Calogero system is still unknown. 
Fourth, can one construct it as a path-ordered exponential of df in a practical way? Fifth, 
what happens if we allow for twisted superfields in the superspace approach? We hope to come 
back to some of these issues in the future. 
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